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The large volume behaviour of the topological susceptibility in 517 (3) gauge theory is investigated on the lattice 
to establish an upper limit on the parity violating terms. 
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1. INTRODUCTION 

A theorem by Vafa and Witten [Q states that 
QCD cannot spontaneously break parity. How- 
ever the validity of this result has been questioned 
by several authors. At finite temperature there 
are parity-odd Lorentz non-invariant operators 
that may have non vanishing expectation values 
in the vacuum It is also not clear whether 
the conclusion of the theorem can be extended to 
purely fermionic operators |2j ■ On the other hand 
the proof followed by the authors of pQ might im- 
plicitly assume the absence of spontaneous sym- 
metry breaking of parity 0]. 

Recently Aguado and Asorey |5I6| have proved 
that parity cannot be broken by the topological 
charge operator Q. The vacuum energy density 
E(Q) is a well defined smooth function for any 
value of 9. This excludes the possibility of a cusp 
at 9 = and implies (Q) = and a volume- 
independent finite value for the topological sus- 
ceptibility x- 

A direct check of this issue on the lattice can 
be useful. The idea is to study the volume 
(in)depcndence of the topological susceptibility 
X- One can parametrize the breaking of parity 
by writing for the topological charge 



Q = ±aV + g 



(1) 



X 



(2) 



where a is a parameter, V = L the space-time 
volume (L is the lattice size). The quantity ±aV 



is the expectation value of Q in each of the two 
(possible) broken parity vacua among which, at 
finite volume, the system can do tunneling (a = 
means no parity violation), g is a fluctuation 
whose distribution satisfies (g) = and (g 2 ) ^ 0. 
Assuming that g is stochastically independent of 
the broken parity vacuum, the topological suscep- 
tibility reads 

V V 
Lattice simulations can give an upper bound to 
\a\ and hence to the possibility of a spontaneous 
breaking of parity. 

2. THE LATTICE PROCEDURE 

Pure Yang-Mills SU{2>) theory was numerically 
simulated by using the standard plaquette action 
on the lattice at /3 = 6.0 on three different vol- 
umes: 16 4 , 32 4 and 48 4 . In all cases a heat- 
bath algorithm was used for updating configura- 
tions and care was taken to decorrelate successive 
measurements in order to render them indepen- 
dent. On the lattice we determined the topolog- 
ical charge with the 1-smeared operator in- 
troduced in [7j (see also [5]). The relationship 
between the lattice bare topological susceptibil- 
ity xl = ((Ql^ )/V an d tne physical one \ IS 
given by 

(3) 



XL = a 4 Z 2 X + M 



1 



2 



where a is the lattice spacing, Z and M are renor- 
malization constants that can be extracted by us- 
ing a nonperturbative method [10111) . In Eq.(3) 
X is the physical topological susceptibility that 
coincides with d 2 E(9)/dd 2 . As it is argued in 
if parity is broken a volume dependence of x is 
expected. 

In order to provide a clear upper bound to the 
parity breaking parameter \a\, we have performed 
high statistics runs for all lattice sizes by using 
the APEmillc computers in Pisa. The results are 
shown in Fig.l. The trend in this figure is consis- 
tent with a volume-independence of xl for large 
enough sizes. However, any definite conclusion 
must be drawn from the physical topological sus- 
ceptibility x- 
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Figure 1. Bare topological susceptibility xl at 
(3 = 6.0 for several lattice volumes. 



If the same values of Z and M are used for all 
volumes, it is clear from Eq.(3) that Fig.l will 
be transformed into an analogous figure that dis- 
plays the same trend for the physical topolog- 
ical susceptibility x- However in perturbation 
theory the renormalization constants present a 
mild dependence on the lattice size L, at least 
for small L. This is a well-known result for 



several 2-dimensional models |12|. Although no 
similar analysis have been carried through for 4- 
dimensional QCD, it is expected that in this the- 
ory the inclusion of the zero-modes will show up 
again as an explicit dependence on L. 

In principle an L dependence is expected also 
for the exact nonperturbative determination of Z 
and M. We are currently calculating Z and M 
following the method of |10lll| . High statistics 
is required to reduce the errors. We are using 
lattice sizes 8 4 , 12 4 and 16 4 . The results will be 
extrapolated by using the functional forms 
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M(L) 
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